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<N 1 Introduction 

S . 

Consider the Schrodinger equation in M 3 with the compactly supported potential q (x) , x G 
M 3 . The problem of the reconstruction of the function q (x) from measurements of the 
solution of that equation on a certain set is called "inverse scattering problem". In this 
paper we prove uniqueness theorems for some 3-d inverse scattering problems in the case 
when only the modulus of the complex valued wave field is measured, while the phase is 
unknown. This is the phaseless case. In the past, phaseless inverse scattering problems were 
studied only in the 1-d case (section 1.2). As to the 3-d inverse scattering problems in the 
frequency domain, it was assumed in all studies so far that both the modulus and the phase 
of the complex valued wave field are measured, see, e.g. [2] for uniqueness results in the case 
of a piecewise analytic potential and (261 for global uniqueness results and reconstruction 
>• . methods. 

Below C s+a are Holder spaces, where s > is an integer and a G (0, 1) . Let fl,Gcl 3 
be two bounded domains, Q C G. For an arbitrary point y G M 3 and for an arbitrary 
number uj G (0, 1) denote B w (y) = {x : \x — y\ < u>} and P w (y) = M. 3 \B LC [y) . For any two 
sets M, N C M 3 let dist (M, V) be the Hausdorff distance between them. Let Gi C M 3 be a 
convex bounded domain with its boundary S G C 1 . Let e G (0, 1) be a number. We assume 
that ficGiCG, dist (S, dG) > 2e and dist (S, <9ft) > 2e. Hence, 

dist(dB E (y),dty > s,VyeS, (1) 
5-i i dist (dB £ (y) , dG) > e, \/y G S. (2) 



Below either m = 2 or m = 4, and we will specify this later. We impose the following 
conditions on the potential q (x) 

q (x) G C m (M 3 ) , q (x) = for x G R 3 \G, (3) 
q(x) > 0. (4) 

As a rule, the minimal smoothness of unknown coefficients is not the first priority of 
proofs of uniqueness theorems of multidimensional coefficient inverse problems, see, e.g. 
[26| [27] and Theorem 4.1 in [29]. Since our proofs require either C 2 or C 4 smoothness of 
solutions of Cauchy problems for some hyperbolic equations, we are not concerned below 
with minimal smoothness assumptions. In particular, the reason of imposing C 4 (rather 
than C 2 ) smoothness conditions in Theorems 3 and 4 is rooted in smoothness requirements 
of uniqueness theorems of [U El EJ [191 1201 [21] , which we use here. 
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1.1 One of main results 



We now formulate one of our four main theorems. Three other theorems are formulated in 
section 2. Let x = (x ,i, 20,2, #0,3) be the source position. Consider the following problem 

A x u + k 2 u — q (x) u = —5 (x — x ) , x G R 3 , (5) 

u (x, xq, k) = O ( -) , \x\ — > 00, (6) 

\\x- Xq\J 

~r ^rd x .u (x, xq, k) — iku (x, xq, k) = o ( -] , \x\ — > 00. (7) 

^—^ x — x \\x — Xq\ J 

Here the radiation conditions ([6]), ([7]) are valid for every fixed source position xq. To establish 
existence and uniqueness of the solution of the problem (J5])-(J7]), we refer to Theorem 6 of 
Chapter 9 of the book [32] as well as to Theorem 3.3 of the paper |31j. As to the smoothness 
of the solution of the problem ([5])-(I7j), we refer to Theorem 6.17 of the book [12J. Thus, 
combining these results, we obtain that for each pair (k, xq) e IxK 3 there exists unique 
solution u (x,xo, k) of the problem ([5]), ([6]), ([7]) such that 

U (x, Xq, k) = Uq (x, X , k) + U s (x, Xq, k) , (8) 

no = ex ^\ x ~ x ^\ Us e cm+i+a {Pu ( Xo) ) , Va,a; G (0, 1) . (9) 

47T \X — Xq\ 

In ([8]), ([9]) Mo (x, xo, k) is the incident spherical wave and u s (x, xo, k) is the scattered wave. 

Inverse Problem 1 (IP1). Let m = 2 in ([3]). Suppose that the function q (x) satisfying 
([3]), (jlj) is unknown for x G Q and known for x G M 3 \fi. Also, assume that the following 
function / x (x, x , k) is known 

fi (x, Xq, k) = \u (x, Xq, k)\ 2 , Vx G S, Vx G B £ (x ) , x ^ x , V/c G (a, b) , 

where (a, b) C M is an arbitrary interval. Determine the function q (x) for i£ll. 
Theorem 1 is one of four main results of this paper. 

Theorem 1. Consider IP 1. Let two potentials qi (x) and q 2 (x) satisfying conditions (T3|] ; 
^ be such that q\ (x) = g 2 [%) = Q (^) for x G M 3 \fi. Let U\ (x,x , k) and u 2 (x,x , k) be 
corresponding solutions of the problem (TJP-{7P satisfying conditions $8$), (TJJ). Assume that 

\ui (x, Xq, k)\ 2 = \u 2 {x, Xq, k)\ 2 , Vx G S, Vx G B £ (x ) , x 7^ x , VA; G (a, b) . (10) 

Then qi (x) = q 2 (x) . 

Corollary 1. Fix two arbitrary points y G S and y G B £ (y ) such that y ^ y . Suppose 
that all conditions of Theorem 1 are in place, except that is replaced with 

\ui(y,yo,k)\ 2 = \u 2 (y,y ,k)\ 2 , \/k G (a, b) . (11) 

Then m (y, y , k) = u 2 (y, y , k) for all k G R. 
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Remark 1. The proof of Corollary 1 can be immediately derived from the proof of 
Theorem 1. Completely analogous corollaries are valid for each of Theorems 2-4 of this 
paper. Their proofs can also be immediately derived from proofs of corresponding theorems. 
We omit formulations of those corollaries for brevity. 

We now outline the main difficulty, which did not allow to prove uniqueness results for 
phaseless 3-d inverse scattering problems so far. As an example we consider IP1. Analogous 
difficulties take place for three other inverse problems formulated in section 2. In IP1 one 
should work with a complex valued function r (k) ,k G R such that its modulus \r (k)\ is 
known for all k G (a,b). The function r (k) admits the analytic continuation from the real 
line R in the half-plane {k G C : Im k > —7} for a certain number 7 > 0. Since \r (k)\ 2 = 
r (k) r (k) , then the function \r (k)\ 2 is analytic for k G R as the function of the real variable 
k. Here r(k) is the complex conjugate of r (k) . Hence, the modulus \r (k)\ is known for all 
k G R. Denote C + = {k G C : Imk > 0} . Proposition 4.2 of |17| implies that if r [k) would 
not have zeros in C+, then this function would be uniquely reconstructed for all k G R from 
the values of |r(fc)| for k G R, also see Lemma 4 in subsection 3.2. However, the main 
difficulty is to properly account for zeros of r (k) in the upper half-plane C+\R. Indeed, let 
zi, z n G C + \R be some of such zeros of r (k) . Consider the function r (k) defined as 



Hence, |r (k)\ = \r (k) \ , V7c G R. Furthermore, the function r (k) is analytic in the half-plane 
{k G C : ImA; > —7} . Therefore, in order to prove uniqueness, one needs to figure out 

how to combine the knowledge of \r (k)\ for k G R with a linkage between the function r (k) 

and the originating differential operator. 

This difficulty was handled in |16j in 1-d, using the fact that the function r (k) depends 

only on one variable k in this case. Unlike [16], the function r (x, xo, k) = u (x, xo, k) depends 

on x,xo,k in the 3-d case. Hence, the above zeros depend now on both x and xo, i.e. 

) . Thus, compared with the 1-d problem, the main difficulty of the 3-d case is 

that it is necessary to figure out how to take into account the dependence of zeros Zj (x, x ) 

from x and Xq. To do this, we essentially use here properties of the solution of the Cauchy 

problem for an associated hyperbolic PDE. 

1.2 Published results 

The phaseless inverse scattering problem is of central importance in some applications, where 
only the amplitude of the scattered signal can be measured. An example is neutron specular 
reflection, see, e.g. [3]- Uniqueness of the phaseless inverse scattering problem in the 1-d case 
was first proved in [16]. Next, the result of [16] was extended to the discontinuous impedance 
case in [23]. Also, see [I] for a relevant result. A survey can be found in [T7]. Uniqueness 
theorem for a 1-d phaseless inverse problem arising in crystallography was proven in [T3]. 
This problem is essentially different from the one considered in [16J. 
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Inverse problems without the phase information are well known in optics, since it is often 
impossible to measure the phase of the optical signal, unlike its amplitude. In optics, such a 
problem is usually formulated as the problem about the recovery of a compactly supported 
complex valued function from the modulus of its Fourier transform. The latter is called the 
"phase retrieval problem" [9]. This problem arises in x-ray crystallography [22], astronomical 
imaging [10J and other subfields of optics [9]. Some numerical methods for this problem can 
be found in, e.g. [H HI [101 EES EQ]. Recently regularization algorithms were developed in the 
1-d case for a similar, the so-called "autocorrelation problem" 0HT]. Uniqueness theorems 
for the phase retrieval problem can be found in [TJ1 [18]. 

In section 2 we formulate three more phaseless inverse scattering problems as well unique- 
ness theorems 2-4 for them. In section 3 we prove Theorem 1. Theorem 2 is proved in section 
4. Finally, Theorems 3 and 4 are proved in section 5. 

2 Other problems and results 

In IP1 the modulus of the total wave field u = uq + u s is known on a certain set. We now 
consider the case when the modulus of the scattered wave is known. 

Inverse Problem 2 (IP2). Let m = 2 in ([3]). Suppose that the function q (x) satisfying 
([3]), (HI) is unknown for x G Q and known for x G M 3 \fi. Also, assume that the following 
function f 2 (x,Xo,k) is known 

f 2 (x,x ,k) = \u a (x,x Q , k)\ 2 , Vx G S,Vx G B £ (x ) ,x ^ x ,Vk G (a, b) . 

Determine the function q (x) for x G fl. 

Theorem 2. Consider IP2. Assume that all conditions of Theorem 1 hold, except that 
(TJ2P is replaced with 

|u S) i (x,x , k)\ 2 = \u s>2 {x,x , k)\ 2 ,Wx G S, Wx G B e (x ),x ^ x ,Wk G (a, b) , (12) 

where u s j = Uj — u ,j = 1,2. In addition, assume that q(x) ^ 0, Va; G S. Then q 1 (x) = 
q 2 (x) . 

Theorems 1 and 2 are formulated only for the over-determined data. Indeed, in both IP1 
and IP2 the number of free variables in the data exceeds the number of free variables in the 
unknown coefficient. The reason of this is that even if the phase is known, still all current 
uniqueness results for 3-d inverse scattering problems in the case when the 6— function is 
the source function are valid only if the data are over-determined ones, see, e.g. [21 [261 121] 
for the frequency domain and §1 of chapter 7 of [21] for an inverse scattering problem in 
the time domain. Suppose now that the function 5 (x — xq) in (jSJ) is replaced with such a 
function p (x) that p (x) ^ in Q. And consider the inverse problem of the reconstruction of 
the potential q (x) from values of the function u (x, k) for all x G S, k G R. Then uniqueness 
theorem for this problem can be proved for the non-overdetermined case. This proof can be 
handled by the method, which was introduced in the originating paper Also, see, e.g. 



4 



[U [HI [20] and sections 1.10, 1.11 of [I] for some follow up works of authors of [5] on this 
method; a survey can be found in [21J. This technique is based on Carleman estimates. 

Consider the function \ (x) G C°° (R 3 ) such that x { x ) — 1 i n G\ an d X ( s ) = f° r x ^ G. 
Let xo £ 5*. For a number a > consider the function 5 CT (x — xo) , 

, X(«) ( \ x ~ x o\ 2 

O a (X — Xq) = C 3 exp 



where the number C > is such that 

5 a (x — xq) dx = 1. 



G 



The function 5 a (x — xo) approximates the function 5 (x — xq) in the distribution sense for 
sufficiently small values of a. The function S a (x — x$) is acceptable in Physics as a proper 
replacement of 5 (x — Xq), since there is no "true" delta-function in the reality. On the other 
hand, the above mentioned method of (5] is applicable to the case when 5 (x — xq) is replaced 
with 5 a (x). Therefore, it seems to be worthy from the Physics standpoint to consider Inverse 
Problems 3,4 below. 

Let in (j3J) m = 4. To apply results, which follow from the method of [5], consider the 
function g (x) such that 

geC 4 (R 3 ) ,g(x)=0m R 3 \G, (13) 
g(x)^0mGx. (14) 

Consider the following problem 

Av + k 2 v - q (x) v = -g (x) ,x G M 3 , (15) 



x\ — > oo, (16) 









, \ x \ 


H 











, \x\ — y oo 







(17) 

The same results of [121 EH [32] as ones in subsection 1.1 guarantee that for each fcefi there 
exists unique solution v (x, k) G C 5+a (M 3 ) , Vet G (0, 1) of the problem (P), dTTJ). 

Inverse Problem 3 (IP3). Let m = 4 in (j3J). Suppose that the function g (x) satisfying 
conditions (jl]) is unknown for x G Q and known for x G M 3 \fi. Assume that the following 
function fs (x, fe) is known 

f 3 (x,k) = \v(x,k)\ 2 ,Vx eS,Vke (a,b) . (18) 

Determine the function q (x) for x G fl 



5 



Theorem 3. Consider IPS. Let the function g (x) satisfies conditions (TO)), JT^j j. Con- 
sider two functions q\ (x) , q 2 {x) satisfying conditions (TJ|) ; and such that q\ (x) = q 2 (x) = 
q(x) for x G R 3 \fi. For j = 1,2 let Vj (x,k) G C 5+a (IR 3 ) be the solution of the problem 
(E5p-|77p with q(x) = qj (x). Assume that 

\vi (x, k) | 2 = \v 2 (x, k) \ 2 , Vx G S, \/k G (a, b) . (19) 

T/ien qi (x) = q 2 (x) . 

We now pose an analog of IP2. Let v (x, k) be the solution of the problem (fTBTl - ffTTI) for 

the case q (x) = 0, 



/" exp(ifc|z -£|) 



Hence, one can interpret the function t> (x, k) as the solution of the problem (fT5|) - (|T7|) for 
case of the background medium. 

Inverse Problem 4 (IP4). Let m = 4 in ([3]). Suppose that the function q (x) satisfying 
(EJ) , (SJ) is unknown for x G Q and known for x G R 3 \f2. Let v s (x, k) = v (x, k) — vq (x, k) . 
Assume that the following function / 4 (x, k) is known 

fi (x, k) = \v s (x, k) | 2 , Vx G 5*, Vfc G (a, 6) . 

Determine the function q (x) for x G fi. 

Theorem 4. Consider IP4- Let all conditions of Theorem 3 hold, except that / TiPj) zs 
replaced with 

\v sA ( Xj fc) | 2 = | Vj>j2 fc) | 2 ? Vx G 5, VA; G (a, 6) , (20) 

where v s j(x,k) = Vj(x,k) — Vq(x, k),j = 1,2. Assume that q(x) ^ 0,Vx G S. Then 
qi (x) = g 2 (#) • 

3 Proof of Theorem 1 

3.1 Functions U and u 

Consider the solution U(x,xo,t) of the following Cauchy problem 

U u = A X U - q (x) U, (x, t) G M 3 x (0, oo) , (21) 

17 (x, 0) = 0, U t (x, 0) = 5 (x - x ) . (22) 
It was shown in §1 of Chapter 7 of the book [24J that the function U has the form 



5(t- 


|x - x |) 


4tt . 


X — X \ 





U (x, x , t) = U \ J± + 17 (x, x , t) , (23) 
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where 

&(w)=- / tmsgjzL^m* , (24) 

D(x,XQ,t) 

D{x,x ,t) = {£: | x -e| + |x -e| <t}, (25) 
U{x,x Q ,t) = for i G (0, [ac — ar [) , (26) 



[/(x,x ,t) = -T^n 1 / q(Z)ds + 0(t-\x-xo\),t->\x-x \ + . (27) 

16n\x — x \ J 

L(x,xo) 

Here L (x, xq) is the interval of the straight line connecting points x, Xq- 

Let $ C M. 3 be an arbitrary bounded domain. Choose the number ujq = ojq (xq, $) G (0, 1) 
so small that $ fl P w (scq) 7^ 0. Let T > max ie 5 \x — Xq\ be an arbitrary number. Denote 

*{$,x ,u,T) = {(x,t) : xe^nP w (x Q ),t G [\x-x \,T]}. 

It was shown in §3 of chapter 2 of the book [28] that the function U can be represented as 
the following series 

00 

U (x, Xq, t) = } j U n (X, Xq, t) , (28) 
n=0 

where 

U (x,x ,t) = 5{t ~ lx ~ X f\ (29) 



U n (x,xo,t) = - / vs ; | 1 ^ de,n>l. (30) 

J 47r|x-£| 

D(x,xo,t) 

Convergence estimates of §3 of chapter 2 of [28] imply that series ([28"]) converges in the 
norm of the space C 2 (\& ($, x , w, T)). Hence, for any fixed source position x G M 3 and for 
|0|=O,1,2 

D^ t U(x,x ,t) GC(f($,%u,T)),Vu6 (0,w (2o,$)),VT>max|2-x |. (31) 

We now refer to some results about the asymptotic behavior of solutions of hyperbolic 
equations as t — > oo. More precisely, we refer to Lemma 6 of Chapter 10 of the book [32] as 
well as to Remark 3 after that lemma. It follows from these results that there exist numbers 
C\ = C\ (q, $, Xq, u) > 0,ci = ci(g, $,xo,w) > depending only on the function q, the 
domain the source position xq and the number u G (0, ojq (xq, $)) such that 

D^ t U(x,x ,t) < C ie - Clt in {(x,t) : x G $ fl P w (xq) , t > |x-x |}, \(3\ = 0,1,2. (32) 
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It follows from (123|) . (!3T|) and (|32|) that one can apply the operator J 7 of the Fourier 
transform with respect to t to functions D^ t U (x, Xq, t) , \/3\ = 0, 1, 2. Let 



7 (U) (x, x Q , k) = U (x, x Q , t) e dt, Vx, x 6l,i^ x , Vk G 



Using again the same results of references [T2], [31], [32] as ones cited in section 1.1, we obtain 
that 

u (x, xq, k) = T (U) (x, xo, k) , Vx, x$ G R 3 , x ^ xo, V/c G R. (33) 

In particular, (1321) and (1331) imply that for each pair x, xo G R 3 such that x 7^ xq there 
exists a number 7 = 7 (x, xo, g) > such that the function u (x, Xq, k) admits the analytic 
continuation with respect to k from the real line R in the half-plane {fceC: ImA; > —7}. 

Suppose that G n P u (x ) ^ 0. Then, using ([23>(j2Tjh (1521 and ([33]), we obtain the 
following asymptotic behavior of functions u (x, x , k) and w s (x, x , fc) for every fixed source 
position xq G 



° 3 and for every fixed value of u G (0, (xo, $)) 
exp (i/c |x — xq|) 



u (x, x , k) 



Air |x — x | 



\k\ — >■ 00, k E C- 



u s (x,x , fc) 



i exp (2A; |x — x |) 
167r |x — xq| /c 



L(x,xo) 



— > OO, Aj G C + , 



(34) 



(35) 



uniformly for x E G C\ (xo) . Therefore, we have proven Lemma 1. 

Lemma 1. The solution U (x,Xo,t) of the problem l[21\) . UM) can be represented in 
the form l[2^) , where the function U(x,xo,t) satisfies conditions (24\)- f3l^) . Furthermore, 
holds, where the function u (x,xo, k) is the unique solution of the problem |3]) ; (0|), (?|) 
satisfying conditions (TJ]) ; (G|). In addition, for every pair of points x,x G R 3 such that 
x ^ xq the function u (x,x , k) admits the analytic continuation with respect to k from the 
real line R in the half-plane {k G C : Imk > — c{\ , where C\ = C\ (q, $,Xo,w) > is the 
number in / fffjgj) . Finally, if G fl P w (xo) 7^ 0, then asymptotic formulas (34\ ) and l[35\) hold 
uniformly for x G G fl P w (xq). 



3.2 Lemmata 2-5 

Lemma 2. Lei Xo G R 3 and x G G, x 7^ Xo &e two arbitrary points. Then the function 
u (x, Xo, k) has at most finite number of zeros in C + . 

Lemma 3. Let xq G R 3 and x EG be two arbitrary points. Assume that 

L(x,x ) 



Then the function u s (x, xq, k) has at most finite number of zeros in C+. 

Lemmata 2 and 3 follow immediately from (|34|) and (|35|) respectively. 

Lemma 4. Let 7 > be a number. Let the function d (k) be analytic in {k G C : Im k > —7} 
and does not have zeros in C+. Assume that 



d(k) 



C_ 

k n 



[l + o (1)] exp (ikL) , \k\ -»■ 00, fc G C H 



where C G C and n ; L G R are some numbers and also n > 0. Tnen £/ie function d(k) 
can be uniquely determined for k G {/c G C : Im > —7} fry £/ie values of \d(k)\ for fceR. 
Furthermore, for fcef 



axgdik) = — lim lim 
7T R->oo £ ->o+ 



k—e 



-R 



I n Id (01 



ln|d(Q| 



72 7T 



+ Lfc-argC + — . (36) 



The right hand side of (1361) can be any of branches of the function arg. However, 
this does not make any difference for us, since we are interested in the function d (k) = 
\d(k) \ exp [i argd (k)] . Since the function d (k) is uniquely determined on the real line, then 
the analyticity of this function implies that it is uniquely determined in {k G C : Im k > —7} . 
Lemma 4 follows immediately from Proposition 4.2 of [IT]. Hence, we omit the proof. 

Lemma 5. Let the function d(k) be analytic for all k G R. Then the function \d(k)\ 
can be uniquely for all k G R by the values of \d (k)\ for k G (a, b). 

This lemma was actually proven in subsection 1.1 for the function r (fc). 



3.3 Proof of Theorem 1 

Choose an arbitrary point Xq G S and an arbitrary point x G B £ (xq) such that x ^ i . 
Denote 

hi (k) = ui (x, x , k) , h 2 (k) = u 2 (x, x , k) . (37) 

It follows from Lemma 1 that we can regard hi (k) and h 2 (k) as analytic functions in the half- 
plane {k G C : ImA; > —61} , where 61 > is a certain number. Next, since \hi (k) \ = \h 2 (k)\ 
for k G (a, 6) , then Lemma 5 implies that 

\h 1 (k)\ = \h 2 (k)\,Vk eR. (38) 

By Lemma 2 each of functions hi (k) , h 2 (k) has at most finite number of zeros in C+. Let 
{ai, a n } C (C + \R) and {61, b m } C (C + \R) be sets of all zeros of functions hi (k) and 
h 2 (k) respectively in the upper half-plane. Here and below each zero is counted as many 
times as its order is. Let {a' 1; a' ri } C R and {b[, ...,&£. 2 } C R be sets of all real zeros of 
functions hi (k) and h 2 (k) respectively. 
First, we prove that 

K ...,<} = {^,...,6^}. (39) 
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Let, for example a[ be the zero of the order m > 1 of the function hi (k) as well as the 
zero of the order mi > of the function h 2 (k) . Then hi (k) = (k — a[) ni hi (k) and hi (k) = 
(k — ai) mi hi (k) , where 

h x {a' x )h 2 {a' x ) (40) 

Using (1381) . we obtain 



(k-a[) ni • /n(jfe) = (ife-ai) mi • M*0 , VA;GR. 



(41) 



Assume, for example that mi < ri\. Dividing both sides of ( 14TT) by \(k — a'i) mi \ and setting 
the limit at k — > a[, we obtain h 2 (a[) = 0. This contradicts to (1401) . Hence, mi = n\. Since 

a'i is an arbitrary element of the set {a' l5 a^} , then ( l39l) follows. 

Let {ci, c r } C R be the set of all real zeros of both functions h x (k) and h 2 (k) . Define 
functions hi (k) , h 2 (k) as 



hi (*) = h (k) ( n * aj 



V7=l 



fc — a,- 



n 



/l Co 



h 2 (k) = h 2 (k)m k hj 



k - &,■ 

3=1 J 



vS=l 



C.o 



(42) 



(43) 



Hence, /ii (k) and /i 2 (A:) are analytic functions in {k G C : Im A; > —9{\ . Also, these functions 
do not have zeros in C + . Furthermore, it follows from ([3*8]) . (|4*2]) and (|4*3]) and /ii (A;) = 

h 2 (k) , VA; G R. Also, using (134|) . we obtain the following asymptotic behavior of both 
functions hi {k) and h 2 {k) 



hi (k) 



exp (ik \x — Xq\) 



|A;| -> 00, A; G C+, j = 1,2. 



4-7T |x — Xo| k r 

Hence, Lemma 4 implies that hi (k) = h 2 (k) , VA; G R. Hence, (1421) and (1431) lead to 

n , in , 7- 

MA;) n = n ^ VA; G 



=1 fc ~ 6 J 



Hence, 



Rewrite (fUJ) as 



m , , n , 



(44) 



hi (k)+hi (k) mt—^.- 1 



1 k — bj 

VJ=1 



^ (A;) + /i 2 (A;) ]"[ 



A; — ay 
A; — cL 



(45) 
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Consider the function w\ (k) , 

k — 0,4 



This function can be rewritten in the form 



n 1 

Wl (k) = Q (k) J] — — 



where Q (k) is a polynomial of the degree less than n. By a partial fraction expansion, A (k) 
can be written in the form 



Wl (k) = £ 



Ci 



where Cj G C are some numbers and Sj,n' > 1 are integers. Direct calculations verify that 
the inverse Fourier transform J 7 " 1 of the function (k — (ij)~ Sj is 

T ~ X ( (fc-g.)* ) = H (t) d ^ exp (-^^ 
with a certain constant Cj G C. Hence, 

n' 

J- 1 (wi (k)) := Ai (t) = H(t)J2 Cjt Si ~ l exp (-iajt) , (46) 

3=1 

where constants Cj = CjCj G C. Similarly, denoting 

k — bj 



7=1 j 



jr-i ( W2 (A;)) : = A 2 (t) =H{t)J2 B P~ X ex P H¥) > ( 47 ) 



we obtain 



where constants B i G C. For j = 1,2 let J 7-1 (fy) = h 5 it) . By ([33]) and ([3j 

/ij (t) = Z7j (x,x ,t) , (48) 

where Uj (x, xq, t) is the solution of the problem ( 12~TT) . f l22l) with g (x) = qj (x) . We now apply 
the operator J 7-1 to both sides of ( )45l) . Using ( 1461) . ( 1471) and the convolution theorem, we 
obtain 

t t 
hi(t)+ h l (t-s)X 2 (s)ds = h 2 (t)+ [h 2 (t-s)X 1 (s)ds. (49) 
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Since Xq G S, x G B £ (xq) , then (JTj), ([2J and (1251) imply that 

D(x,x ,t) c ,Vt G (|x-x | (50) 

Since q x (x) = q 2 (x) for x G M 3 \ft, then ([251), (125} , (p]) and (j5U]) imply that 

Ui (x, Xq, t) = U 2 (x, Xq, t) = U (x, Xq, t) , Wt G (\x — X \ , s) . (51) 

By (gED and (EU) 

hi (t) = h 2 (t) = h{t) = U (x, x , t),Vte(\x-x \,e). (52) 
Hence, PHI) implies that 



t 



J h(t-s)\(s)ds = 0, Vt G (|a; -%| ,e) , (53) 
o 

where 

A (t) = Ax (t) - A 2 (t) . (54) 
Using (1231), (126), (12TI), (l3Tj) and (152)1, we obtain 

An \x — Xq\ 



p(t) = 0,tE(0,\x-x \), (56) 
1 

p (t) = ■ 

t^\x-x \+ 167r|x — Xq\ 



lim p (t) = — r / q (f ) ds, (57) 



L(x,xo) 

•i2 



peC 2 (t>\x-x \). (58) 

Introduce a new variable r t, where r = t — \x — x \ . Then (1531) and (I55l) - (l58l) lead to the 
following integral equation of the Volterra type with the continuous kernel p (r + \y — y \ — s) 
and with respect to the function A (r) 

T 

A (r) + An \x — xq\ J p (r + \x — x \ — s) A (s) = 0, Vr G (0, £ — |x — x \) . (59) 
o 

Hence, A (r) = for all r G (0, £ — \x — x \) . On the other hand, (j4"B"l) . (|47|) and (154)) imply 
that A (t) is analytic function of the real variable t > 0. Hence, A (£) = 0, Vt > 0. This implies 
that {ai, a n } = {bi, b m } . Thus, (1371) and (14*41) lead to u\ (x, Xq, k) = u 2 (x, xq, k) , V7c G 
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Since Xo G S and x G B £ (xq) are two arbitrary points such that x 7^ Xq, then we have 
established that 



Mi (x, x , k) = u 2 (x, Xq, k) , Vx G S, Vx G -B e (x ) , x 7^ x , VA; G M. 



(60) 



Consider an arbitrary point y G 5. Since qi (x) = g 2 (x) = g (x) in M 3 \f2, then, using (J5]), 
we obtain 

A x Uj + fc 2 Mj — q (x) Mj = — 5 (x — ?/ ) j x ^ M 3 \fi, j = 1, 2. 

Hence, ( 1601 and the well known result about uniqueness of the continuation problem for 
elliptic equations (see, e.g. §1 of chapter 4 of [21]) imply that u\ (x, yo, k) = u 2 (x, yo, k) , Vx G 
M 3 \fi, VA; G R. Hence, 



Mi (x, Xo, A;) = M2 (x, xo, A;) , Vx, xo G S", x 7^ xo, VA; G M. 
Using (|g|l. (15}. (13"41) and (1331). we obtain for j = 1, 2 



lim 

fc— ><x> 



II ; 



4ik ( — — 1 ) (x, x , k) 

M 



g.,- (£) ds, Vx, x G M 3 , x 7^ x . 



(61) 



(62) 



L(x,xq) 



Hence, (JSP) and lead to 



j (<7i - 92) (0 rfs = 0, Vx, x G 5. 



(63) 



L(x,x ) 



Finally, (163]) and the classical uniqueness theorem for the Radon transform implies that 
<?i (x) = q 2 (x) ■ □ 

The idea of using (I63p for the proof of the uniqueness of an inverse scattering problem in 
the time domain can be found in §1 of chapter 7 of [23] 



4 Proof of Theorem 2 

Consider an arbitrary point xo G S. Since q\ (x) = q 2 (x) = q (x) for x G !R 3 \f2 and 
q (x) ytz o, Vx G S, then ([1]) and (T5]) imply that one can choose e > so small that 

J qi(0ds= J q 2 (0ds= J q (£) ds ± 0,\/y G B £ (x ) . (64) 

L{y,x ) L(y,x ) L(y,x ) 

Choose an arbitrary point x G B £ (xo) such that x 7^ Xo- Denote h s j (k) = u s j (x,xo, k) . It 
follows from (I64|) and Lemmata 1, 3, 4 and 5 that we can apply to functions h s j (k) the same 
procedure as the one described in section 3.3. For brevity we keep notations (|46|) . (|47|) and 
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(1531). For j = 1,2 let /i s>i (i) = J 7 " 1 (/i Sji ) . Let Z7(x,ar ,t) be the function defined in ([23]) . 
(I24|) . Then, using (!55|) . ( |56|) and ( |58|) . we obtain analogously with (152]) 

(t) = £7 (x,x ,t) = p(t) ,Vt € (|x - x | ,e) ,j = 1,2. (65) 
Hence, using (1651) . we obtain similarly with (1591) 



p (r — s + \x — Xq\) A (s) ds = 0, Vr G (0, e — \x — x |) , (66) 



where r = t — \x — xq \ . Differentiating both sides of (|66|) with respect to r and using ([5 
([58]) and (J64J), we obtain 



A (r) — m (x, x ) J p' (r + \x — x \ — s) X (s) ds = 0, Vr G (0, e — |x — x |) , 

Q 

m(x,x ) = I 16?r I / <?(0 rfs 

\L(x,x ) ) 

Hence, similarly with the proof of Theorem 1, we conclude that A (t) = 0, Vt > 0. This leads 
to 

u St i (x, x , k) = u Si2 (x, x , k) ,\/k G R. (67) 

Next, since Uj (x, xo, k) = u (x, xo, k) + u S) j (x, xq, k) , j = 1, 2 and since again xo <E S and 
x G B £ (x ) are two arbitrary points such that x ^ x , then (IBTj) implies (1601) . The rest of 
the proof is the same as the proof of Theorem 1 after (1601 . □ 

5 Proofs of Theorems 3 and 4 

5.1 Functions V and t> 

Consider the following Cauchy problem 

V tt = AV-q (x) V, {x, t) G R 3 x (0, oo) , (68) 

V (x, 0) = 0, V t (x, 0)=g (x) . (69) 
Corollary 4.2 of Chapter 4 of the book [23] implies that there exists unique solution V G 



H 2 (R 3 x (0,T)) , VT > of the problem (EHD, (J69J). By the Kirchhoff formula this function 
V (x, t) is the solution of the following integral equation 

^ CM) = 1=7 / 9(.£)dSz- / — r- „ (70) 



4?rt 
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Construct functions V n (x, t) as 



V (x,t) = ± j g(Z)dSt, 



\x-t\=t 



V n {x,t) 



q(OV n -i &t~ 






x-£\ 



n = 1,2,..., 



(71) 



(72) 



\x-£\<t 



The above mentioned technique of §3 of chapter 2 of the book [28] implies that the function 
V (x, t) can be represented as 



V(x,t) = Y i V n (x,t) 



(73) 



n=0 



and this series converges in the norm of the space C 4 ($ x [0,T]) for any bounded domain 
Id 3 and for any number T > 0. Hence, 



V G C 4 ($ x [0,T]) 



(74) 



Using again Lemma 6 in Chapter 10 of the book [32] as well as Remark 3 after that lemma, 
we obtain that for any bounded domain $ C M. 3 there exist constants C% = C 2 (q,g,<&) > 
0, C2 = C2 (q, g, $) > depending only on functions q, g and the domain $ such that 



D% t V(x,t) < C 2 e~ C2 \ Vx G $,Vt > 0; = 0,1,..., 4. 



(75) 



Hence, using again Theorem 6 of Chapter 9 of [32], Theorem 3.3 of [31] and Theorem 6.17 
of [T2] . we obtain that the Fourier transform of the function V (x,t) is the unique solution 
v (z, fe) G C 5+a (M 3 ) , Va G (0, 1) of the problem (TT5l)-(fT7l). i.e. 



v{x, k) = F{V) ,Vx G M 3 ,VA; G M. 



(76) 



Furthermore, it follows from (1751) and (1751) that for every point x G G the function u (2, fc) 
admits the analytic continuation with respect to k from the real line K in the half-plane 
{k G C : Im k > -c 2 (q, g, G)} . 

Hence, using the integration by parts in integral (!76|) of the Fourier transform as well 
as (|69|) . (!74|) and (!75|) . we obtain the following asymptotic behavior of the function v (x, k) , 
uniformly for x G G 



v {x, k) = ^ 



-g(x)+0[- 



\k\ — > 00, k G C+. 



Let 
Then 



K(x,t) = V(a;,*) - VoOM). 
w s (x, fc) = F{V S ) . 



(77) 

(78) 
(79) 
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Next, using (168]) and ( 169|) . we obtain 

d 2 t V s = AV S - q (x) (V s + V ) , (x, t) G R 3 x (0, oo) , (80) 

V s (x, 0) = d t V s {x, 0) = 0. (81) 
Hence, ([55]), (ISO]) and (1ST]) imply that 

<9[K (x, 0) = for r = 0, 1, 2 and d 3 V s (x, 0) = - (g#) (x) . (82) 

Hence, using ( 182]) and the integration by parts in the right hand side of (1791) . we obtain the 
following asymptotic behavior of the function v s (x, k) , uniformly for x G G 

v s (x, k) = ^[- (qg) (x) + o (1)] , \k\ -» oo, k G C+. (83) 
Thus, we have proven Lemma 6. 

Lemma 6. There exists unique solution V (x, t) of the problem / TffPj) such that (7~4\ ) 

is valid for every bounded domain $ C M. 3 and for every T > 0. Estimate ( [73] ) is valid for 
this function V (x, t) . In addition, the Fourier transform [76] ) of the function V (x, t) is the 
unique solution v(x,k) G C 5+a (R 3 ) , Vet G (0,1) of the problem (T73])-([77]). Also, for every 
x G $ the function v (x, k) admits the analytic continuation with respect to k from the real 
line R in the half-plane {k e C : Imk > — c 2 } . Finally, asymptotic formulas (77] ) and [8^1 
hold uniformly for x G G. 

Lemma 7 follows immediately from Lemma 6, ( TT4]) . ( 177]) and (183]) . 

Lemma 7. For ei>en/ point x G Gi £/iere exists at most finite number of zeros of the 
function v (x, k) in C + . Next, assume that there exists a point x' G G\ such that q (x') ^ 0. 
Then the function v s (x',k) has at most finite number of zeros in C + . 

5.2 Proof of Theorem 3 

Consider an arbitrary point x G S and denote similarly with (137]) hj (k) = Vj (x, k) ,j = 1, 2. 
By Lemma 6 there exists a number 62 > such that each of functions h± (k), hi (k) admits 
the analytic continuation from the real line M in the half-plane {k G C : Imk > —82} ■ By 
Lemma 7 each function hi (k) , hi {k) has at most finite number of zeros in C+. Hence, the 
asymptotic behavior (]77j) enables us to apply the technique of section 3.3. For j = 1, 2 let 
hj (t) = J 7-1 (hj) . Let the function Vj (x,t) satisfying (JTJJ be the solution of the problem 
(168]) . (169]) with q (x) := qj (x). Then (176]) implies that 

hj(t) = Vj(x,t),j = 1,2. (84) 
Since x e S, then it follows from and ([70])- ([73]) that 

Vi (x, t) = V 2 (x, t) = V (x, t) , Vt G (0, e) . (85) 
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We briefly note that another way of establishing (!85|) is via the energy estimate. Using ([8] 
and (185|) . we obtain 

h x (t) = h 2 (t) = V {x, t) , Vt G (0, e) . (86) 

Using arguments, which are completely analogous with those of section 3.3, keeping the same 
notations (jlSj) . (jUJ), (j5Ij) and using (IBS"]) , we obtain the following analog of (166]) 

t 

U (x, t — t) A (r) g?t = 0, Vt G (0, e) . (87) 



o 

Differentiating both sides of ( IHTj) twice, using initial conditions ( 169]) as well as ( 174]) . we obtain 
the following Volterra integral equation of the second kind with respect to the function c (t) 

t 

A(t) + -Lr / U tt (x,t - r) A (r) rfr = 0, Vt G (0,e) . (88) 
o 

Hence, analogously with subsection 3.3, we conclude that A (t) = 0, Vt > and 

vt {x, k) = v 2 {x, k),Vke E. (89) 

Since x G S is an arbitrary point, then (|89j) implies that 

U! (x, jfe) = w 2 {x, k) , Vx G 5, VA; G E. (90) 

Since Vj (x, t) = J 7-1 (uy) , j = 1, 2, then ( 190]) leads to 

ViOM) |5 00 =^ 2 (x,t) \ Soo := v (x,t), (91) 

where S'oo = S x (0, oo) . Let V (x, t) be any of two functions V\ (x, t) , V2 (x, t) . Since 

gi (a?) = q 2 (x) = q (x) for x G M 3 \fi, (92) 

then (168]) . (169]) and f]9T]) imply that the function V (x, t) satisfies the following conditions 

V tt = AV-q (x) V", (x, t) G (M 3 \Gi) x (0, 00) , (93) 

V (x, 0) = 0, d t V (x, 0)=g( x ),xe R 3 \Gi, (94) 
^(^*) \ Sac =v(x,t). (95) 

Recall that both functions Vj G (M 3 x (0,T)) , VT > 0. On the other hand, the stan- 
dard energy estimate tells us that the problem (T93])-(f95]) has at most one solution V G 
H 1 ((M 3 \G!) x (0, T)) , VT > 0. Hence, V 1 (x, t) = V 2 (x, t) = V (x, t) for (x, t) G (R 3 \(?i) x 
(0,oo). 
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Denote if (x,t) = d u V (x,t) {s^, where v is the unit normal vector at S pointing outside 
of the domain G\. For j = 1,2 consider functions Wj (x, t) = d t Vj (x, t) . Then 

d$Wj = AWj - qj (x) Wj, (x, t)eG x x (0, oo) , (96) 

Wj (x, 0)=g (x) , d t Wj (x, 0) = 0, (97) 

Wj{x,t) \ Soc = Vt (x,t), (98) 

d u Wj{x,t) \ Soo =<Pt(x,t). (99) 

In addition, (I7i|) implies that 

Wj G C 3 (H x [0, T]) , VT > 0. (100) 

Finally, using Theorem 4.7 of [19], we obtain that relations fl96|) -f ll00|) imply that qi (x) = 
q2 (x) in Q. To finish the proof, we refer to ( 1921) . □ 

Remark 2. Let the number T > diam(G\) /2. Then one can replace in ( 196]) the time 
interval (0, oo) with (0, T) and also replace in ( 1981) . ( 199]) with St = S x (0, T) . These new 
conditions (j55]l . (15811 and (1991 together with (ETJ) and (TTOOj) still imply that q 1 (x) = q 2 (x) in 
Q, see either Theorem 1.10.5.2 of [I], or Theorem 3.2.2 of [20J, or Theorem 3.2 of |21j . 



5.3 Proof of Theorem 4 

Consider again an arbitrary point x G S. For j = 1,2 let hj (k) = v S j (x,k) and hj (t) = 
T v (hj) . Using ( 17X1) and ( 1781) . we obtain similarly with (185!) 

£i (*) = ^2 (*) = K (x, i) , V* G (0, e) . (101) 

We again apply arguments, which are completely analogous to those of section 3.3. Also, 
we keep the same notations (PIS"]) . (I4"7j) and (154]) . Hence, we obtain from (jlOip the following 
analog of (!H7|) 

y K (x, t-r) \ (r) dr = 0, Vt G (0, e) . (102) 
o 

Differentiate (I102p four times and use (1741) and (j82p . We obtain the following analog of the 
Volterra integral equation (1881) 

t 

X{t)- r 4r^ I d?V s (x,t-T)\(r)dT = 0, Vte(0,e). 

(<?#) 0) y 

o 

Hence, similarly with subsection 3.3 we conclude that A (t) = 0, Vt > 0, which implies that 
(x, k) = v S)2 (x, k) ,Wk e JR.. Since x 6 5 is an arbitrary point, then 

t» Sjl (x, jfe) = v s , 2 (x, k) , Vx G S, VA; G R. (103) 

Since Uj (x, A;) = t> s j (x, fc) + vq (x, fc) , then (11031) implies that (1901) is valid again. The rest 
of the proof is the same as the proof of Theorem 3 after (190]) . □ 



18 



Acknowledgment 

This research was supported by US Army Research Laboratory and US Army Research 
Office grant W911NF-11-1-0399. 

References 

[1] T. Aktosun and P.E. Sacks, Inverse problem on the line without phase information, 
Inverse Problems, 14, 211-224, 1998. 

[2] Ju. M. Berezansky, The uniqueness theorem in the inverse problem of spectral analysis 
for the Schrodinger equation, American Mathematical Society Translations, Series 2, 
Volume 35, 137-235, 1964. 

[3] N.F. Berk and C.F. Majkrzak, Statistical analysis of phase-inversion neutron specular 
reflectivity, Langmmr, 25, 4132-4144, 2009. 

[4] L. Beilina and M.V. Klibanov, Approximate Global Convergence and Adaptivity for 
Coefficient Inverse Problems, Springer, New York, 2012. 

[5] A.L. Bukhgeim and M.V. Klibanov, Uniqueness in the large of a class of multidimen- 
sional inverse problems, Soviet Math. Doklady, 17, 244-247, 1981. 

[6] A.L. Bukhgeim, Introduction in the Theory of Inverse Problems, VSP, Utrecht, 2000. 

[7] Z. Dai and P.K. Lamm, Local regularization for the nonlinear autoconvolution problem, 
SI AM J. Numerical Analysis, 46, 832-868, 2008. 

[8] D. Dobson, Phase reconstruction via nonlinear least squares, Inverse Problems, 8, 541- 
557, 1992. 

[9] J.R. Fienup, Phase retrieval algorithms: a personal tour [invited], Applied Optics, 52, 
45-56, 2013. 

[10] J.R. Fienup, J.C. Maron, T.J. Schulz and J.H. Seldin, Hubble space telescope charac- 
terization by using phase retrieval algorithms, Applied Optics, 32, 1747-1768, 1993. 

[11] D. Gerth, B. Hoffman, S. Birkholz, S. Koke and G. Steinmeyer, Regularization of an 
autoconvolution problem in ultrashort laser pulse characterization, Inverse Problems in 
Science and Engineering, to appear. 

[12] D. Gilbarg and N.S. Trudinger, Elliptic Partial Differential Equations of Second Order, 
Springer, New York, 1984. 

[13] N.E. Hurt, Phase Retrieval and Zero Crossings: Mathematical Methods in Image Re- 
construction, Kluwer Academic, Dodrecht, 2002. 



19 



[14] M.V. Klibanov, Determination of a function with compact support from the absolute 
value of its Fourier transform, and an inverse scattering problem, Differential Equations, 
22, 1232-1240, 1987. 

[15] M.V. Klibanov, Uniqueness of the determination of distortions of a crystal lattice by 
the X-ray diffraction method in a continuous dynamical model, Differential Equations, 
25, 520-527, 1989. 

[16] M.V. Klibanov and P.E. Sacks, Phaseless inverse scattering and the phase problem in 
optics, J. Math. Physics, 33, 3813-3821, 1992. 

[17] M.V. Klibanov, P.E. Sacks and A.V. Tikhonravov, The phase retrieval problem. Topical 
Review. Inverse Problems, 11, 1-28, 1995. 

[18] M.V. Klibanov, On the recovery of a 2-D function from the modulus of its Fourier 
transform, J. Mathematical Analysis and Applications, 323, 818-843, 2006. 

[19] M. V. Klibanov, Inverse problems and Carleman estimates, Inverse Problems, 8, 575- 
596, 1992. 

[20] M.V. Klibanov and A. Timonov, Carleman Estimates for Coefficient Inverse Problems 
and Numerical Applications, VSP, Utrecht, 2004. 

[21] M.V. Klibanov, Carleman estimates for global uniqueness, stability and numerical meth- 
ods for inverse problems, J. Inverse and Ill-Posed Problems, 21, issue 2, 2013, to appear; 
available online at arxiv 1210.1780vl. 

[22] M.F.C. Ladd and R.A. Palmer, Structure Determination by X-Ray Crystallography, 
Plenum Press, New York, 1993. 

[23] O.A. Ladyzhenskaya, Boundary Value Problems of Mathematical Physics, Springer, New 
York, 1985. 

[24] M.M. Lavrentiev, V.G. Romanov and S.P Shishatskii, Ill-Posed Problems of Mathemat- 
ical Physics and Analysis, AMS, Providence, R.I., 1986. 

[25] Z.T. Nazarchuk, R.O. Hryniv and A.T. Synyavskyy, Reconstruction of the impedance 
Schrodinger equation from the modulus of the reflection coefficients, Wave Motion, 49, 
719-736, 2012. 

[26] R.G. Novikov, Multidimensional inverse spectral problem for the equation — Aip + 
(v(x) — Eu(x))ijj = 0, Functional Analysis and Its Applications, 22, 263-272, 1988. 

[27] R.G. Novikov, The inverse scattering problem on a fixed energy level for the two- 
dimensional Schrodinger operator, J. Functional Analysis, 103, 409-463, 1992. 



20 



[28] V.G. Romanov, Integral Geometry and Inverse Problems for Hyperbolic Equations, 
Springer - Verlag, Berlin, 1974. 

[29] V.G. Romanov, Inverse Problems of Mathematical Physics, VNU Science Press, Utrecht, 
1986. 

[30] B. Sixou, V. Davidoiu, M. Langer and F. Peyrin, Absorption and phase retrieval with 
Tikhonov and joint sparsity regularizations, Inverse Problems and Imaging, 7, 267-282, 
2013. 

[31] B.R. Vainberg, Principles of radiation, limiting absorption and limiting amplitude in 
the general theory of partial differential equations, Russian Math. Surveys, 21, 115-193, 
1966. 

[32] B.R. Vainberg, Asymptotic Methods in Equations of Mathematical Physics, Gordon and 
Breach Science Publishers, New York, 1989. 

Michael Victor Klibanov 

University of North Carolina at Charlotte 

Department of Mathematics and Statistics 

Charlotte, NC 28223 

U.S.A. 

mklibanv@uncc.edu 



21 



